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ABSTRACT 

We present optical, radio and X-ray data that confirm a new supernova remnant (SNR) 
in the Large Magellanic Cloud (LMC) discovered using our deep Ha imagery. Optically, 
the new SNR has a somewhat filamentary morphology and a diameter of 56 x 64 
arcsec (13.5 x 15.5pc at the 49.9kpc distance of the LMC). Spectroscopic follow-up 
of multiple regions show high [Sn]/Ha emission-line ratios ranging from 0.66 ±0.02 
to 0.93 ±0.01, all of which are typical of an SNR. We found radio counterparts for 
this object using our new Australia Telescope Compact Array (ATCA) 6 cm pointed 
observations as well as a number of available radio surveys at 8 640 MHz, 4850 MHz, 
1 377 MHz and 843 MHz. With these combined data we provide a spectral index a ~ - 
0.5 between 843 and 8 640 MHz. Both spectral line analysis and the magnetic field 
strength, ranging from 124 — 184 /uG, suggest a dynamical age between ^2,200 and 
~4,700 yrs. The SNR has a previously catalogued X-ray counterpart listed as HP 483 
in the ROSAT Position Sensitive Proportional Counter (PSPC) catalogue. 

Key words: Supernova remnants, Magellanic Clouds, ISM: individual objects: MC¬ 
SNR J0512-6707; RP1577, Surveys 


1 INTRODUCTION 

SNRs represent a prime source of heavy elements which are 
synthesised in their high mass progenitors. They therefore 
play a major role in the enrichment of their host galaxies. 
Supernova (SN) explosions produce shocks, powerful blast 
waves and a radiation front strong enough for the resulting 
SNRs to shape and heat the interstellar medium (ISM), im¬ 
pacting on the ionisation balance of the entire area in which 
they reside. SNRs usually produce radiation across the whole 
electromagnetic spectrum with detections in the optical, in¬ 
frared (Lakicevic et al. 2015), radio (due to non-thermal ra¬ 
diation e.g. Stupar et al. 2007), X-ray (detected by satellites 
e.g. Aschenbach 2002), the region of soft gamma-rays (due 
to the presence of residual young neutron stars) and even in 
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the ultraviolet (Blair et al. 2006). Multi-wavelength studies 
of SNRs provide information on several crucial aspects of 
stellar evolution and the physical behaviour of shocked ma¬ 
terial. They have a major influence on local magnetic fields 
whilst the shock waves may trigger new star formation and 
accelerate cosmic rays. 

Within SNRs there is a correlation between Ha and 
radio-continuum emission even though the detailed fine 
structures do not overlap very often (see Cram, Green & 
Bock, 1998). If the SN blasts into a hot, low-density ISM, 
Ha emission may be faint but if the ISM is cool and dense, 
Ha radiation is usually detected (with an intensity depen¬ 
dent on the density and temperature at the interface) reveal¬ 
ing the shape and density of the shock front in the form of 
narrow filaments, elongated arc structures, bubble-like fea¬ 
tures and knots. These same structures can usually be seen 
in the forbidden line emission from [Sn] 6717/6731 and to 
a more varying degree in [N il] 6548/6584, [O i] 6300/6364, 
[O III] 4959/5007 and the [O il] doublets close to 3727 A and 
7325 A. Following Fesen, Blair & Kirshner (1985), a ratio 
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Figure 1 . The median stacked and combined Ha and broad band 
red optical image of RP1577 where Ha is shown as red and the 
short red is shown as white. 

of [S n] to Ha > 0.5 implies a shock front indicative of an 
SNR as does strong [O i] and [O n] lines plus the supporting 
presence of [N n] (Sabbadin et al. 1977). In addition, the de¬ 
tection of Ha, rather than [O ill], leading the radiation front 
would support the rejection of any Wolf-Rayet (WR) shell 
scenario (Gruendl et al. 2000). 

At a distance of ~49.9kpc (Reid & Parker 2010 and 
references therein) the LMC is an ideal laboratory to study 
the different varieties of SNRs in great detail across an entire 
galaxy. The LMC is close enough that SNRs are resolved suf¬ 
ficiently for multiwavelength studies to be conducted. Since 
SNRs yield key information on the death rate of high mass 
stars and the number of white dwarfs exploding as Type 
la SN, an improved inventory of LMC SNRs, complete as 
possible, is required to put constraints on these values and 
improve our understanding of the final stages of stellar evo¬ 
lution. The LMC SNR sample is currently estimated at 59 
remnant^] 

The deep, stacked UKST Ha survey used to search 
for planetary nebulae in the LMC (Reid & Parker, 2006a, 
2006b, 2013), also revealed several candidate SNRs, the first 
of which we report here, thanks to the optical depth of the 
Ha survey to R eq uiv mag 22, and the method of image anal¬ 
ysis. Spectroscopic followup of these objects (originally iden¬ 
tified as SNR candidates in 2004) using 2dF on the Anglo- 
Australian Telescope (AAT) and the 1.9-m telescope at the 
South African Astronomical Observatory (SAAO) have per¬ 
mitted us to confirm identifications and conduct an optical 
spectroscopic and multiwavelength analysis. 


1 estimated by the authors using 31 SNRs from Chandra cat¬ 
alog (hea-www.cfa.harvard.edu), 54 from Badenes et al. (2010) 
and recent discoveries including those from Maggi et al. (2014), 
Bozzetto et al. (2014a,b), Warth et al. (2014) and Kavanagh et 
al. (2015) 


The rest of the paper is structured as follows. In sec¬ 
tion [2] we present the optical observational data including 
images and spectra. In subsection |2. 1 1 we derive an estimate 
for the expansion velocity and age of the SNR and in |2.2| 
estimate the radial velocity from key emission lines. In sec¬ 
tion [3] we provide supporting radio data, the results of in¬ 
frared imaging and previous X-ray observations. Section [4] 
is our conclusion where we summarize the combined multi¬ 
wavelength data used to achieve our classification. 


2 OPTICAL OBSERVATIONS AND 

DISCOVERY 

The new SNR candidate, given the Reid & Parker refer¬ 
ence name RP1577 and presented here for the first time, 
was discovered using deep-stack Ha and short red images 
of the LMC, SMC and surrounding regions (Reid & Parker 
2006a). An Ha interference filter of exceptional specification 
and quality (Parker & Bland-Hawthorn 1998) was used to 
cover a substantial fraction of the UKST’s large field. Its 
central wavelength (6590 A) and bandpass (70 A FWHM) 
work effectively in the UKST’s fast f/2.48 converging beam. 
With a peak transmission of >85 per cent, the Ha filter 
is also partially sensitive to [Nii] 6583A. Overlapping fields 
with non-standard 4-degree centres enabled full contiguous 
coverage in Ha despite the circular aperture. High resolu¬ 
tion, panchromatic Tech-Pan film with peak sensitivity at 
Ha was used as the photographic detector (see Parker et al. 
2005 for more details). 

Over a period of three years, from 1997, a series of 12 
repeated narrow-band Ha and 6 matching broad-band ‘SR’ 
(Short Red) exposures of the central LMC field were also 
taken (see Reid & Parker 2006a). The purpose was that these 
images should be stacked to produce a single Ha and short 
red image with an overall depth increase of 1 magnitude over 
a single image frame. The twelve highest quality and well- 
matched UK Schmidt Telescope 2-hour Ha exposures and 
six 15-minute equivalent SR-band exposures were selected. 
From these exposures, deep, homogeneous, narrow-band Ha 
and matching broad-band SR maps of the entire central 25 
deg 2 square of the LMC were constructed. 

The ‘SuperCOSMOS’ plate-measuring machine at the 
Royal Observatory Edinburgh (Hambly et al. 2001) scanned, 
co-added and pixel matched these exposures creating 10/mi 
(0.67 arcsec) pixel data which goes 1.35 and 1 magnitudes 
deeper than individual exposures, achieving the full canon¬ 
ical Poissonian depth gain, e.g. Bland-Hawthorn, Shopbell 
& Malin (1993). This gives a depth ~21.5 for the SR images 
and Requiv ~22 for Ha (4.5 x 10 _17 er<?s cm~ 2 s _1 A -1 ) 
which is at least 1-magnitude deeper than the best wide-field 
narrow-band LMC images previously available. An accurate 
world co-ordinate system was applied to yield sub-arcsec as¬ 
trometry, essential for success of the spectroscopic follow-up 
observations. 

The resulting Ha and short red (SR) maps were origi¬ 
nally overlayed with false colours in order to detect plane¬ 
tary nebulae (Reid & Parker, 2006a,b,2013) but also proved 
an excellent tool for uncovering a large number of hot and 
cool emission-line stars (Reid & Parker, 2012). In addition to 
these discrete objects, the maps were able to reveal the fine 
morphological characteristics of SNRs, which appear as com- 
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Table 1. Observing Logs for the follow-up of emission object and SNR candidate RP1577 in the LMC. The first column gives the 
telescope name, second gives the instrument name, third gives the observation date, the fourth gives the grating name, the fifth gives 
the dispersion in angstroms, the sixth gives the central wavelength, the seventh gives the wavelength coverage in angstroms, the eighth 
gives the exposure time in seconds, the ninth gives the number of exposures taken. 


Telesc. 

Inst. 

Date 

Grating 

Dispenser 

Dispersion 

A/pixel 

Central 

A (A) 

Coverage 

A (A) 

T exp 

S 

Neorp 

AAT 

2dF 

16 Dec-04 

300B 

4.3 

5852 

3660 - 8000 

1200 

3 

AAT 

2dF 

18 Dec-04 

1200R 

1.1 

6793 

6220 - 7340 

1200 

2 

1.9-m 

ccd 

7 March-14 

300 

5 

5800 

3850 - 7738 

1800 

2 


plex interwoven streams, elongated arcs, bubbles and narrow 
filaments in Ha. The newly identified candidate is shown in 
Fig. [l] where Ha is coloured red and the short red continuum 
image appears white due to both blue and green colours be¬ 
ing assigned to the same image. Although these structures 
can be closely correlated with radio emission, they usually 
become quite fragmented in old remnants and may extend 
beyond, within or along the boundaries of the radio remnant 
(see Stupar & Parker 2009). In addition to the Ha/SR com¬ 
bined images, we used the Magellanic Cloud Emission-Line 
Survey (MCELS; Smith et al. 1998) to directly compare the 
strength of the three narrow bands of [O ill] 5007 (blue), Ha 
(red) and [Sn] 6716+6731 (green). When the colour maps 
are evenly combined the SNR candidate becomes evident 
by the strong green ([S il] line) colour. 

RP1577 has a central position of (J2000) RA 05 h 12 m 
28 s DEC -67° 07' 20” determined using SuperCOSMOS as¬ 
trometry (Hambly et al 2001) and a measured angular extent 
from the Ha imagery of 56 x 64 arcsec. This equates to a 
physical extent of 13.5 x 15.5 pc at the LMC distance. With 
physical diameters ranging from 0 to 120 pc (Badenes et al. 
2010), the cumulative distribution of SNR diameters in the 
LMC have been found to be roughly linear (eg. Chu & Ken- 
nicutt 1988; Bandiera & Petruk 2010) up to a cutoff at a 
physical radius of r cu ^ ~30pc (Badenes et al. 2010). This 
means that the 15.5 pc size of RP1577 places it among the 
smaller 15 per cent of LMC SNRs, joining six other SNRs 
in the 10-20 pc range. 

To eliminate the possibility that RP1577 was the prod¬ 
uct of a WR shell we checked the available databases for the 
presence of any WR star or candidate inside (or in the vicin¬ 
ity of) the shell but none were found within a three aremin 
radius. We then compared the expanded radius of Ha and 
[O ill] emission since, according to theoretical models, the 
outermost [O ill] usually leads the Ha emission in WR shells 
(Gruendl et al. 2000). This occurs due to a decrease in the 
sensitivity of Ha at temperatures greater than 10 4 K while 
[O ill] emission is unaffected by the higher temperatures. As 
a result, the morphologies of SNRs and WR shells should be 
somewhat different (see for example Cox 1972, Gruendl et 
al. 2000). As these models depend on the interaction of the 
shells with the interstellar medium, they remain as a use¬ 
ful diagnostic whether or not they are used to identify WR 
shells. Using the MCELS survey we overlayed the [O ill] and 
Ha maps in order to check whether [O ill] was leading at any 
point around RP1577. Despite the low 2.5 arcsec/pixel res¬ 
olution of the MCELS survey, a pixel by pixel comparison 
along the shock front showed that Ha was in fact leading 
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Figure 2. The (Ha/SR) quotient image of RP1577 showing cir¬ 
cular loops to the north of parallel, stratified filaments which 
outline the exterior of the expanding shell. The position of the 
2dF spectroscopic observation is shown by the yellow circle close 
to the center of the object. Horizontal lines mark the positions of 
long-slit observations. 



Figure 3. The 2dF spectroscopic observation of RP1577 indi¬ 
cated by the yellow circle in Figure [ 2 ] The Ha/[Sn] ratio is 
0.46zb0.02. The extinction of c(H/3) = 0.45 is based on the Balmer 
decrement. 
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Table 2. Observed line intensities and flux ratios for four different spectroscopically observed positions across three areas of RP1577. 
Line measurements were made using the splot script in the iraf package. Errors for individual flux line measurements are a constant 
±1.9xl0 -17 for 2dF 300B, ±1.2xl0~ 17 for 2dF 1200R, ±1.5xl0 -17 for AAOmega, ±2xl0 -17 for 1.9-m measurements. Column details 
are provided in the text. 


Obs. 

(1) 

RA 

( J2000) 

(2) 

DEC 
(J2000) 

(3) 

H/3 

(4) 

Flux (erg < 
[O in]5007 
(5) 

=ms-l) 

Ha 

(6) 

[ N n ]6583 
(7) 

[Nil]/ 

Ha 

(8) 

[S ii] / 

Ha 

(9) 

[S II] ratio 
6716/6731 
(10) 

(=m 3 ) 

(11) 

T e 

(K) 

(12) 

2dF-300B 

05 12 28.4 

-67 07 29 

2.22E-15 

3.84E-15 

6.95E-15 

1.49E-15 

0.32±0.02 

0.46±0.02 

1.35±0.03 

70±9 

35,000± 10,000 

2dF-1200R 

05 12 28.4 

-67 07 29 

- 

- 

- 

- 

0.25±0.02 

0.66±0.02 

1.44±0.03 

<26 

- 

1.9-m 1 

05 12 31.0 

-67 06 59 

1.81E-14 

2.20E-14 

6.40E-14 

1.16E-14 

0.24±0.01 

0.78±0.01 

1.40±0.02 

26.5±8 

- 

1.9-m 2 

05 12 26.0 

-67 07 05 

3.2E-15 

7.17E-15 

1.17E-14 

2.32E-15 

0.30±0.01 

0.93±0.01 

1.44±0.02 

<26 

- 




Figure 4. The 2dF 1200R medium resolution spectroscopic ob¬ 
servation of RP1577 indicated by the yellow circle in Figure [2] 
The higher resolution of this spectrum shows the weak [O i]6300 
& 6363 lines. The [Sll]/Ha ratio is 0.66+0.02. 


Figure 5. The long-slit, spectroscopic observation of RP1577 at 
position (J2000) RA 05 h 12 m 31 s and DEC -67° 06' 59". The 
[S 11 ]/lift observed ratio at this position is 0.78=b0.01. The extinc¬ 
tion of c{ H/3) = 0.28 is based on the Balmer decrement. 


the [O ill] emission, clearly supporting the rejection of any 
WR shell scenario. 

In Table [l] we summarise details regarding the spectro¬ 
scopic follow-up observations. A five night observing run on 
the AAT using 2dF (Lewis et al. 2002) was undertaken in 
December 2004 to spectroscopically confirm LMC emission 
candidates. The large corrector lens incorporates an atmo¬ 
spheric dispersion compensator, which is essential for wide 
wavelength coverage using small diameter fibres. Specifically 
targeted long-slit spectra were obtained using the 1.9-m tele¬ 
scope at the SAAO in March, 2014. 

Optical properties for the new SNR in the LMC are 
given in Table [5] Column (1) indicates the telescope and 
observing method used at each position. The AAT observa¬ 
tions were conducted using the 2dF multi-fibre spectroscopy 
system where a each fibre covered 2 arcsecs diameter on the 
sky. Column (2) gives the (J2000) position of the optical cen¬ 
ter of the observed position in right ascension. Column (3) 
likewise, gives the (J2000) declination of the same position. 
Each observed position was determined using the Ha map 
and checked against SuperCOSMOS astrometry. A compar¬ 
ison of the resulting spectral signature was intended to pro¬ 
vide a good overall analysis of the ratio of [Sn] to Ha and 
[N n] across the nebula, all measured with gaussian line fit¬ 
ting from within the iraf splot script. Column (4) gives the 
measured H/3 flux at the position of each slit or fibre. Column 
(5) gives the measured flux of the [O ill] 5007 emission line 
at the observed position. Flux calibration of the 2dF line in¬ 
tensity was conducted using the method outlined in Reid & 



Figure 6. The long-slit, spectroscopic observation of RP1577 at 
position (J2000) RA 05 h 12 m 26 s and DEC -67° 07' 05". The 
[SllJ/Ha observed ratio of 0.93+0.01 at this position is the high¬ 
est ratio of all our spectroscopic observations of RP1577. The ex¬ 
tinction of c(H/3) = 0.33 at this position is based on the Balmer 
decrement. 


Parker (2010). For the purpose of diagnostics the Ha flux is 
provided in column (6), followed by the [N n] flux in column 
(7). Column (8) then gives the [N il] 6583+6548/Ha ratio 
including the error based on line measurements. It should 
be noted that the [N n] 6548 line is quantum mechanically 
fixed at a 1-3 intensity ratio with the [N n] 6583 line. Col- 
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umn (9) provides the [Sll]/Ha ratio and line measurement 
errors for estimating the excitation of the shock-front. Col¬ 
umn (10) shows the ratio of [S n] 6716/6731 used for esti¬ 
mating the electron density. The electron density shown in 
column (11) is derived from the [S n] 6716/6731 ratio derived 
using the temden script found in the stsdas package from 
iraf. Once the [Sn] ratio reaches 1.4, representing a den¬ 
sity of 26 cm 3 the low density limit using this method has 
been reached. Column (13) gives the electron temperature, 
derived from the [O ill] 5007+4959/(0 il] 4363 emission lines. 
The faint [O il] 4363 line was only detectable for the purpose 
of measurement at one position. 

The optical Ha/H/3 ratio at each observed position was 
used to determine the extinction constant c(H/J) (i.e., the 
logarithmic extinction at H/l) at that position. These hydro¬ 
gen transitions are the strongest and easiest to accurately 
measure in the nebula spectrum and are fortunately located 
very close to the other main emission lines used for diag¬ 
nostics. The observed Ha/H/3 ratio, when compared to the 
recombination value of 2.86 (Aller 1984), gives a logarithmic 
extinction at H/3 of: 

c(H/3) = (log(Ha/R/3) - log (2.86))/0.34 (1) 

This estimation is based on the relationship between 
observed and intrinsic intensities: 

/ 0 i, s (Ha) _ hatfRaf G (H/9)[/(Ha)-/(H/3)] / 9 s 

I obs ( H/3) Iint(m ’ K 1 

where [/ (Ha)-/ (H/3)] = -0.34 from the standard interstellar 
extinction curve. The resulting extinction values c(H/3) are 
given in the caption below each spectrum in Figs. [3][5] and [6] 

The [Sn] 6717/6731 ratio of 1.35+0.07 from the 2dF 
300B exposures results in an electron density value near 
70 ±10 cm 3 . In Fig. [5] we show the (Ha/SR) quotient image 
which effectively removes stars from the area. The yellow 
circle in the image marks the position observed with 2dF. 
This is approximately 10 arcsec south of the optical centre as 
indicated by the elliptical shape of the diffuse Ha emission. 
Spectroscopic results from this position are shown in Fig. [3] 
for the 300B exposure and Fig. [4] for the 1200R exposure. 

Both 2dF spectra show characteristics and ratios that 
are quite typical of many other SNRs we have observed 
within the LMC. The [Nll]/Ha and [S ll]6716/6731 line ra¬ 
tios from table [2] are in full agreement with those of other 
LMC SNRs (Payne et al. 2008). At the position observed 
by 2dF the [Sll]/Ha ratio of 0.46+0.02 from the 300B grat¬ 
ing is commonly expected for an SNR (Fesen et al. 1985). 
Higher [Sll]/Ha ratios of 0.78 and 0.93 were obtained from 
the 1.9-m spectra in the northern regions of the object. In 
these positions our image data show dense filaments of emis¬ 
sion which can be created in an SNR by post-shock plasma. 
Just below the centre of the SNR, the shocks have a differ¬ 
ent effect, increasing the electron temperature and ionisa¬ 
tion structure of the recombination region, as also seen in 
the more prominent [Ol] 6300 lines in Figs. [ 3 ] and [4] At the 
adiabatic stage in the evolution of the SNR we would expect 
the photoionisation fields to be rather weak close to the cen¬ 
tre where photons are re-ionising gas which is seeping back 
into the shock. This probably explains why the [0 1 ] 6300 
line is higher in the central region but lower in the outer 
shells where the [Sll]/Ha ratios appear to be higher. The 
[Ol]6300 lines in the raw 1.9-m long-slit spectra are wide 


near the continuum therefore their low level in the reduced 
spectra (Figs. |5] and |6| may also in part be due to sky sub¬ 
traction where the combined nebulae and telluric [0 1 ] 6300 
lines span the width of the slit. Morphologically, although 
the overall shape is somewhat elliptical, the interior detail 
of the SNR shows a series of elongated filaments, which be¬ 
come more dense and finally curved in the form of loops to 
the north. With such strong spectroscopic confirmation, we 
now assign RP1577 the name MCSNR J0512—6707 which 
includes the basic coordinates as recommended by the IAU. 

2.1 Expansion velocity and age 

To estimate the expansion velocity and the dynamical age, 
we first assume that the SNR is currently in the Sedov- 
Taylor phase (Sedov, 1959) where the swept-up mass is large 
compared to the ejected mass, the shock velocity is decel¬ 
erating and the evolution is approximately adiabatic. The 
dynamical age is then t = C ■ R/ V s , with C as a constant 
which depends on the expanding velocity of shock wave af¬ 
ter the SN explosion. This method, often used for shell-type 
remnants in the adiabatic stage would require C to be a 
value of 2/5. In this scenario MCSNR J0512—6707 is still 
in the Sedov-Taylor phase, which is typical for SNRs up to 
~30,000yrs of age and diameters up to ~30pc (Badenes et 
al. 2010). As such the radii should expand as 

r ~ El /5 p~ 1/5 t 2/5 (3) 

where Eo is the kinetic energy of the original SN, p is the 
ambient gas density and t is time. The shock velocity decel¬ 
erates according to 

V. = ^ ~ E^p-^r 3 ' 5 (4) 

Our adopted ambient initial gas density is very diffi¬ 
cult to estimate given the location of the SNR very close to 
a dense Hu region and surrounded by dense dust emission 
(see section 3.3). Estimates may be complicated by inhomo¬ 
geneous regions of nebula with varying optical thickness or 
optically thick cloudlets of dust. For a close approximation 
to the average initial H density, we adopt the 1 cm' 3 am¬ 
bient gas density found for use with LMC SNRs by Berhui- 
jsen (1987). This value may be compared to the initial H 
density value of 10cm -3 for Hu regions in the LMC as 
found by Pellegrini et al (2011) which takes into account 
gas-phase abundances and the gas-to dust ratio averaged 
across the densest regions. To estimate an error in the am¬ 
bient gas density we allow for a rise of up to 1.2 cm -3 due 
to dense dust in the immediate area. This was estimated by 
measuring diffuse emission in the general 30 arcmin radius 
surrounding MCSNR J0512—6707 using our 2dF fibres. The 
[S 11 ] 6716,6731 lines provide electron densities (n e ) which are 
proportional to the square root of the emission-line volume 
emissivity. These measurements provide a clue as to whether 
the ambient environment of MCSNR J0512—6707 is propor¬ 
tionally more dense than average. The electron density n e of 
ambient Hu is between 55cm -3 and the low density limit 
of = 26.0 cm -3 which is slightly higher than the standard 
LMC measurement of faint, diffuse emission (< 26.0cm -3 ). 
Since we find that diffuse Hu regions have a low density 
n e = ~600cm -3 , and a H density of 7 cm -3 (see above), 
an n e ~ 38 cm -3 on the same scale would be equivalent to a 
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Table 3. Heliocentric radial velocities and systematic errors for 
key emission lines at the observed positions as shown in Fig. [2] 


Obs. 

[O in] 

5006.84A 

km s -1 

Ha 

6562.80A 

km s -1 

[No] 

6583.4lA 

km s -1 

[So] 

6716.47A 

kms- 1 

error 

km s -1 

2dF-i200R 

_ 

377 

391 

396 

±12 

Posn. 1 

332 

341 

376 

365 

±28 

Posn. 2 

312 

390 

407 

411 

±28 


density of ~1.2 cm^ 3 . For the lower error estimate we adopt 
a value of 0.2 cm -3 since a value of 0.3 cm -3 was found by 
plotting the SNR on the correlation provided by Berezhko 
& Volk (2004). 

The upper error for our adopted energy is estimated 
from the radio observations for this remnant (see section |3| 
and based on a surface brightness of 1.57 x 10 -20 W n? 
Hz -1 sr -1 at a frequency of 1 GHz, placing the object at a 
position on the E - D diagram by Berezhko & Volk (2004) 
equivalent to an energy (Eg) of ~ 0.25 x 10 51 erg. We adopt 
this estimated value as it represents an upper estimate for 
such an SNR in the LMC. For our lower error we adopt 
a value of 9.5xl0 51 erg, which is just below the canonical 
10 51 erg value for SNRs. The radius of 7.75 pc is found 
from the maximum diameter given above. This method pro¬ 
vides an estimated dynamical age of the remnant between 
~2,000yrs and ~5,000yrs. The method assumes that the 
swept-up mass will become proportionally larger than the 
ejected mass and the velocity will decelerate accordingly. 
Using the same equations, the expansion velocity may lie 
between ~ 600 and 1400 km s -1 . 


2.2 Radial velocity 

The heliocentric radial velocities at each observed position of 
the SNR are shown in Table[3] Emission lines were measured 
using the IRAF splot and EMSAO tasks against standard 
rest velocities. Column 1 gives the name of the observa¬ 
tion, which corresponds to the positions marked in Fig. [2] 
Columns 2 to 5 then present the heliocentric radial velocities 
in km s -1 for the [O ill] 5007, Ha, [N il] 6583 and [Sn] 6716 
emission lines. The estimated velocity error, given in column 
6 is a standard systematic error for all velocities measured 
on a particular observation. It was determined according to 
the size of the pixels, related to the resolution of the in¬ 
dividual spectrum. The rest velocities for each element are 
provided in the column headers. 

As measured through optical imagery in section 2, 
above, the [O ill] emission is less affected by temperatures 
greater than 10 4 K and so, can be seen to move at a slower 
velocity than Ha which leads the shock front. This effect can 
be measured, even though expansion at positions 1 & 2 are 
directionally near to the plane of the sky. Our velocity mea¬ 
surements show that this is true of the [O ill] line at both 
positions observed by the SAAO 1.9-m telescope. Velocities 
on the forbidden lines of [N n] and [S n] are in every case 
higher than those found for Ha. This is to be expected for 
[S n] since the bulk abundance of this element is produced 
in the SN explosion. 



5 h 13"30 s 13“00 s 12 m 30 s 12 m 00 s 


Right Ascension (J2000) 

Figure 7. Radio contours from PMN at 4850 MHz super¬ 
imposed on the 843 MHz SUMSS greyscale map of SNR 
MCSNR J0512—6707. Contours are spaced at 0.01 mjybeam -1 . 
The flux density varies from -0.66 mjy to 3.16 mjy across the 
area shown. 

3 RADIO, INFRARED AND X-RAY 

OBSERVATIONS 

3.1 Radio observations 

Radio detection has traditionally been the main method 
for identification and confirmation of SNRs because of their 
non-thermal spectra, radio shell structures and polarization 
(e.g., Bozzetto et al. 2014b). Radio observations also allevi¬ 
ate problems associated with optical detection where dust 
often obscures optical detail. With the development of new 
instrumentation, the sensitivity of existing radio telescopes 
is constantly being enhanced, allowing us to improve our 
knowledge of SNRs within the LMC. For this project we 
have used a combination of dedicated radio observations 
complimented by freely available radio survey data. 

The radio counterpart to this optically detected SNR 
was originally found using data from the Parkes-MIT-NRAO 
(PMN; beam size = 300") radio survey at 4 850 MHz (Con¬ 
don, Griffith, Wright, 1993) and the Sydney University Mo- 
longlo Sky Survey (SUMSS; beam size = 45") at 843 MHz 
(Bock, Large, Sadler, 1999). Although a radio source is de¬ 
tected at the position of MCSNR J0512—6707in both of 
these surveys, it is more extensive in the PMN due to 
the lower resolution of the survey (see Fig. [7|. The con¬ 
tours therefore do not represent the true size. A detection 
of MCSNR J0512—6707 was made in Filipovic et al. (1995; 
listed as LMC B0512-6710) where a radio continuum anal¬ 
ysis of the LMC was performed at various frequencies using 
the Parkes radio telescope. The SUMMS 843 MHz measure¬ 
ment was omitted from spectral index calculations as it was 
superseded by the 843 MHz Molonglo Sythesis Telescope 
(MOST) measurement, due to the higher sensitivity of the 
MOST image and its greater UV-coverage. 

The PMN and SUMSS radio archival data, though help¬ 
ful in confirmation was of insufficient quality to reveal the 
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5 h 13 In 00 5 1 2 m 45 5 12 m 30 s 12 m 15 s 12 m 00 5 

Right Ascension (J3000) 



Right Ascension (J2GQ0) 


Figure 8. MCSNR J0512—6707 with contours from the Dickel et 
al. (2010) 8640MHz mosaic image (beam size = 22 X 22arc.sec). 
The local rms is er = 0.5 mjy beam -1 and contours are 3, 4, 5, 6, 7 
and 8 a. The radio contours are superimposed on the continuum- 
subtracted MCELS map of the immediate area where Ha is red, 
[O III] is blue and [S II] is green. 


structure and so the latest radio-continuum data and our 
own pointed observations were used to provide details and 
additional frequencies to demonstrate the non-thermal na¬ 
ture of MCSNR J0512—6707. Radio-continuum data used 
in this project includes a 36 cm (843 MHz) MOST mo¬ 
saic image (as described in Mills et al. 19841, a 20 cm 
(1377 MHz) mosaic image from |Hughes et ~ah (20071, and 
6 cm (4800 MHz) and 3 cm (8 640 MHz) mosaic images 
from Dickel et al. (2005, 2010) which are shown as contours 
against the continuum-subtracted MCELS three-colour map 
(see section 2) in Fig. [8] 

These mosaic data were complemented by pointed 
observations of MCSNR J0512—6707 from our Australian 
Telescope Compact Array (ATCA) observations taken on 
2011 November 15 and 16, using the Compact Array Broad¬ 
band Backend (CABB) receiver which has a bandwidth of 
2 GHz. The compact array configuration EW367 was chosen 
in an effort to obtain as much emission from the shorter spac- 
ings at 6 and 3 cm (v = 5 500 and 9 000 MHz, respectively) 
as possible. The resulting 6 cm contours are overlaid on the 
MCELS Ha map of MCSNR J0512-6707 in Fig. [9] Total 
integration time over the t wo days for this source amounted 
to ~64min. The MlRlAxj^ (Sault, Teuben & Wrightjil995) 
and karmaP] ( |Goo ch 1995) software packages were used for 
reduction and analysis. Images were formed using miriad 
multi-frequency synthesis ([Sault & Wieringa| l9l)4) and nat¬ 
ural weighting. They were deconvolved with primary beam 
correction applied. 

To estimate the radio spectrum of 


Figure 9. Contours using the entire 2 GHz bandwidth centered at 
5 500 MHz are overlaid on the Ha image of MCSNR J0512—6707. 
The contours are spaced at ~ 0.01 mjy beam -1 commencing at 
0.3 mjy beam -1 . 

MCSNR J0512—6707, integrated flux density measure¬ 
ments were taken from the four mosaics at frequencies 
of 843, 1377, 4 800, and 8 640 MHz, in addition to our 
observations at 5 500 and 9 000 MHz. However, as in¬ 
terferometers such as the ATCA are an arrangement of 
individual antennae, they lack the short spacings that are 
responsible for the large-scale emission (e.g., of objects such 
as SNRs), resulting in a loss in flux density which becomes 
more pronounced at higher frequencies. For this reason, 
the 9 000 MHz integrated flux density measurements were 
discarded. The 5 500 MHz observations were split into 
four 512 MHz bands (4 732, 5 244, 5 756, and 6 268 MHz), 
and integrated flux density measurements were taken per 
band with fitted gaussians using the miriad task IMFIT. 
The results of these measurements and those taken from 
the mosaics are listed in Table [4] and shown in Fig. |10[ 
Errors in these measurements predominately arose from 
uncertainties in defining the ‘edge’ of the remnant. However, 
we estimate errors to be less than 10 per cent. The rms 
noise estimates were taken using the miriad task SIGEST. 
The resulting spectral index between 843 and 8 640 MHz 

of a =-0.52 ± 0.04 is consistent with the typical SNR 

spectral index value of-0.5. 


3.2 Magnetic field strength 

We use t he equipartition formula as given bv |Arbutina et| 
al. (2012 and the corresponding ‘calculatorrj) to estimate 
the magnetic field strength of this SNR. However, further to 
being unable to accurately measure the extent of the rem¬ 
nant, the width of the shell was impossible to measure from 
a source which is somewhat unresolved in the radio, and 


2 http://www.atnf.csiro.au/computing/software/miriad/ 
1 http://www.atnf.csiro.au/computing/software/karma/ 


4 The calculator is available at 

http://poincare.matf.bg.ac.rs/~arbo/eqp/ 
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Table 4. Integrated Flux Densities of MCSNR J0512—6707 used in this work. 


V 

(MHz) 

A 

(cm) 

rms 

(mjy) 

Flux density 
(mjy) 

Beam Size 
(arcsec) 

Image 

Ref. 

843 

36 

0.8 

115.5 

46.4 x 43.0 

MOST 

1377 

20 

1.1 

86.9 

40.0 x 40.0 

Hughes et al. (2007) 

4 732 

6 

0.2 

49.6 

44.0 x 28.4 

this work 

4 800 

6 

0.7 

53.7 

35.0 x 35.0 

Dickel et al. (2010) 

5 244 

6 

0.2 

47.1 

40.1 x 25.9 

this work 

5 756 

6 

0.2 

44.3 

36.6 x 23.7 

this work 

6 268 

6 

0.2 

40.7 

33.9 x 21.9 

this work 

8 640 

3 

0.7 

29.6 

22.0 x 22.0 

Dickel et al. (2010) 


Table 5. Magnetic field strength and minimum energy for varied 
filling factors (ff). 



Figure 10. Radio-continuum spectral index of 
MCSNR J0512—6707 between 843 and 8 640 MHz. The blue 
markers represent flux density error margins of 10 per cent. 


consequently, the filling factor was unable to be estimated. 
Therefore, Tablc[5]shows the average equipartition field over 
the whole shell of MCSNR J0512—6707 in addition to the 
estimated minimum energies based on various filling factors 
ranging from 1 (a completely filled SNR) to 0.25 (a mostly 
empty shell-like remnant). The range of values for the mag¬ 
netic field, from 124 — 184 pG, are relatively high for an 
SNR, and may be indicative of a younger remnant. This 
can be compared with other LMC SNRs using the same 
method, e.g., those listed in Table [6] Such values may sug¬ 
gest a remnant of a few thousand years, which is equal to or 
less than the age estimate of ~2,280 to ~4,690 yrs derived 
in section but still in-keeping with the remnant’s rela¬ 
tively small size which ranks it as the 6th smallest currently 
known in the LMC (Badenes et al. 2010; Bozzetto et al. in 
prep). 


3.3 Dust 

SNRs heat surrounding dust while mass maps indicate that 
they also destroy 3.7 jlj'® Mq °f dust per SNR through the 
process of sputtering (Lakicevic et al. 2015). Type-II SN 


jj 

B 

(pG) 

J-'mxn 

(ergs) 

1.00 

124 

6.75 X 10 49 

0.75 

134 

5.96 x 10 49 

0.50 

151 

5.01 x 10 49 

0.25 

184 

3.71 x 10 49 


Table 6. Magnetic field strength for a sample of LMC SNRs with 
their associated ages. 


SNR Age AAge B Reference 

(yr) (yr) (mG) 


0509-6731 

400* 

120 

168 

Bozzetto 


2014b 

0519-6902 

600* 

200 

171 

Bozzetto 


2012 a 

0533-7202 

17 500 

3 500 

45 

Bozzetto 


2013 

0508-6902 

22 500 

2 500 

28 

Bozzetto 


2014a 

0529-6653 

37500 

12 500 

48 

Bozzetto 


2012 b 


* These age estimates are taken from |Rest et al.| |2005} 


are able to produce dust on the order of a solar mass per 
explosion (Deneault et al. 2003; Nozawa et al. 2003). The 
average lifetime for interstellar dust under exposed to a SNe 
are estimated at 2 l 4 ' 9 years (Lakicevic et al. 2015) indi¬ 
cating that the amount of dust that survives may be ex¬ 
tremely low. Since SNRs have such a major influence on 
the local dust environment and effect a galaxy’s evolution 
(Lakicevic et al. 2015) it is worth examining the environ¬ 
ment of MCSNR J0512—6707 using available MIR imaging. 
The MIR maps used in this analysis are from “Surveying 
the Agents of Galaxy Evolution” (SAGE) conducted by the 
Spitzer Space Telescope (Meixner et al. 2006). The 8 pm 
data is from the Infrared Array Camera (IRAC: Fazio et al. 
2004) and the 24 pm band is from the Multiband Imaging 
Photometer for Spitzer (MIPS; Rieke et al. 2004). The pixel 
sizes are 0.6 " for all IRAC wavelengths and 2.49 " for 24 pm. 

An image from the combined Ha (red), 8pm (blue) and 
24pm (green) maps (Fig. Ill shows dense dust emission to 
the south which is clearly unrelated to the non-thermal SNR. 
Since the strong dust content at 8 pm dominates the other 
three shorter IRAC wave-bands (see Reid, 2014) they have 
not been shown in Fig. EH The image shows a dust filament 
that appears to emanate from the strong emission to the 
south, proceed northwards, along the eastern border of the 
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5 h 12 m 50 s 12 m 40 s 13 m 30 s 12 m 20 s 12 m 10 s 


Right Ascension (J2000) 

Figure 11. The Ho map of MCSNR J0512—6707 and it’s sur¬ 
rounding area (shown in red) is overlaid with the 8//rri SAGE 
image (blue) and the 24/mi MIPS image (green). 

SNR, and then move across its northern boundary to the 
west. It therefore appears that the SNR is encased in dust, 
however, close examination from our line of sight shows a 
faint, diffuse cloud of 24 pun dust across the face of the SNR. 

Although the 24 pirn emission is spread in patches across 
the wider area, the contribution from 8 ptm, shown in blue, 
is more pronounced in the densest areas of 24 ptm dust. The 
24 pirn dust dominates the region immediately surrounding 
the SNR, however, patches of 8 pirn emission are detectable 
to the east and north. Depending on our viewing angle, it is 
possible that the expanding SNR has impacted and cleared 
the 8 ptm dust towards the SW and along the bright, outer 
shell at J2000 RA 05 h 12 m 30 s DEC -67° 07' 55" in the 
North/North-East. A dense region of 8 pim dust can be de¬ 
tected south of the southern perimeter of the SNR which 
is also immediately south of what appears to be the radio 
center (see Figs. 11 

3.4 X-ray 

The presence of X-ray emission is typical of most SNRs and 
is thought to arise from locations where shocks are con¬ 
verting explosive energy into cosmic ray energy. The ra¬ 
diative signatures of GeV and TeV particles indicate ra¬ 
dio to X-ray synchrotron emission from relativisitc elec¬ 
trons as well as 7 -ray emission from accelerated electrons 
and ions (Vink, 2012). The X-ray source corresponding to 
MCSNR J0512—6707 may be found in Haberl & Pietsch 
(1999) ROSAT PSPC catalogue of X-ray sources in the LMC 
region, where it is listed as HP 483. The central coordinates 
of the source are J2000 RA 05^ 12 m 28.0 s DEC -67° 07' 
27'', which is to within ~4 arcsec of the centre of the rem¬ 
nant. The X-ray source extent from the central position is 
given as 8.5±0.2 arcsec. The count rate (PSPC 0.1-2.4keV) 
is given as 3.33x lCV ' 5 ±8.7x 1CV 4 and the hardness ratios 


are given as follows: HI is 1.0 ±0.38 and H2 is 0.04 ±0.23. 
Further X-ray observations will be provided in Kavanagh et 
al. A&A, submitted. 


4 CONCLUSION 

We have discovered a new SNR in the LMC. This classifi¬ 
cation is based on combined evidence from optical imaging 
and spectroscopy, radio detections with a non-thermal spec¬ 
tral index and an X-ray counterpart. Both optical image 
shown in Fig. [l] clearly depict the shell structure and typi¬ 
cal SNR-like morphology. The most distinguishing features 
are a series of bright loops across the northern boundary of 
the object. An [Sll]/Ha ratio of 0.93 at one position across 
the outer shell suggests a typically strong shock. The pres¬ 
ence and strength of the [0 11 ] 3727, 3729 doublet and the 
[O ill] 4959, 5007 lines are typical for SNRs. We note that the 
nitrogen [Nil]/Ha ratios are low compared to those found 
in the spectra of Galactic SNRs, but this can be attributed 
to the fact that this element is less abundant in the LMC 
compared to the Galaxy (see Dopita, 1979). Confusion with 
Wolf-Rayet ejecta can be ruled out as no WR star was de¬ 
tected in the vicinity and the MCELS survey showed Ha, 
not [Oiii], leading the emission. Assuming that the SNR is 
in the adibatic stage, the Sedov-Taylor method was used 
to find the dynamical age of the remnant which may be 
between ~2,000yrs and ~5,000yrs. The expansion velocity 
will therefore be between ~ 600 km s -1 and 1400 km s -1 . 

There is a very strong radio source at the optical cen¬ 
ter of the SNR. Radio-continuum data used in this project 
includes 36, 20, 6 and 3 cm mosaic images. These mo¬ 
saic data were complemented by pointed observations of 
MCSNR J0512—6707 from our ATCA observations. The 
results of these measurements and those taken from the 
mosaics are listed in Table [4] and shown in Fig. [TO] The 
resulting spectral index between 843 and 8 640 MHz of 
a = —0.52±0.04 is consistent with the typical SNR spectral 
index value of —0.5. Values derived for the magnetic field, 
from 124— 184 /rG, are relatively high for an SNR, and may 
suggest an age of only a few thousand years. 

This SNR is surrounded by an area filled with dense and 
often compact dust, much of which may have been expelled 
outwards or cleared by the expanding shell of the SNR. Any 
contribution the SNR has made towards the local dust con¬ 
tent is impossible to judge due to the large region of dense 
dust to the immediate south. An X-ray detection at the cen¬ 
ter of the SNR was made through a previous survey (Haberl 
& Pietsch, 1999). 
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